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ABSTRACT 

We construct models for strongly-magnetized neutron star atmospheres composed of 
mid-Z elements (carbon, oxygen and neon) with magnetic fields B = lO'^-lO'^ G and effec- 
tive temperatures Tg^ = (1 - 5) x 10^ K; this is done by first addressing the physics relevant to 
strongly-magnetized plasmas and calculating the equation of state and polarization-dependent 
opacities. We then obtain the atmosphere structure and spectrum by solving the radiative trans- 
fer equations in hydrostatic and radiative equilibrium. In contrast to hydrogen opacities at 
the relevant temperatures, mid-Z element opacities are dominated by numerous bound-bound 
and bound-free transitions. Consequently, temperature profiles are closer to grey profiles, and 
photosphere densities are lower than in the hydrogen case. Mid-Z element atmosphere spectra 
are significantly softer than hydrogen atmosphere spectra and show numerous absorption lines 
and edges. The atmosphere spectra depend strongly on surface composition and magnetic field 
but weakly on surface gravity. Absorption lines are primarily broadened by motional Stark ef- 
fects and the (unknown) surface magnetic field distribution. When magnetic field variation 
is not severe, substructure in broad absorption features can be resolved by (phase-resolved) 
CCD spectroscopy from Chandra and XMM-Newton. Given the multiple absorption features 
seen in several isolated neutron stars, it is possible to determine the surface composition, 
magnetic field, temperature, and gravitational redshift with existing X-ray data; we present 
qualitative comparisons between our model spectra and the neutron stars 1E1207. 4-5209 and 
RX J1605. 3+3249. Future high-resolution X-ray missions such as Constellation-X mea- 
sure the gravitational redshift with high accuracy by resolving narrow absorption features; 
when combined with radius measurements, it will be possible to uniquely determine the mass 
and radius of isolated neutron stars. 
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1 INTRODUCTION 

Recent observations by Chandra and XMM-Newton detected spec- 
tral features from isolated neutron stars (INSs). A single or multiple 
absorption features have been found from nearby radio-quiet NSs 
I Haberl et alj|2003l: Ivan Kerkwiik et al.ll2004 iHaberl et al.ll2004bl : 



IZane et alllOOllSchwope et al.i2005l : lHaberll2006l) ^ and two broad 
absorptio n features vyere dete c ted from 1E1207 . 4-5209 (hereafter 
1E12 07; ISanwal et al.l |2002| : iMereghetti et"ai] |2002| : iMori et al.l 
l2005l) . In fact, 1E1207 is unique among these INSs since one of the 
features appears above 1 keV, while the other INSs have absorption 
features at E ^ 0.2-0.7 keV. Identification of the observed spectral 
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features will not only allow a measurement of the surface compo- 
sition and magnetic field strength B but also constrain the nuclear 
equation of state via a measurement of the gravitational redshift Zg, 
where Zg is given by (1 + z^) = (1 - 2GM/Rc^y^'^ and M and R are 
the NS mass and radius, respectively. 

The surface composition is usually assumed to be hydro- 
gen because the fast sedimentation time (due to the strong sur- 
face gravity) causes the lightest elements to rise to the top 
jAlcock & Illarionov|[T980l) and only a tiny amount of hydrogen is 
required to produce an optically thick atmosphere Upmani 198^). 
Hydrogen atmospheres in strong magnetic fields have been studied 
in great detail by addressing th e complicated prope r ties of strongly- 
magn e tized dense plasmas jPavlov et ai] 1 19951: Ikai & Salpetd 
1 19971 : IPotekhin et all 1 19991 : IPotekhin & Chabried |2003| , l2004t) . 
Partially-ionized hydrogen atmosphere models show spectral fea- 
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tures, such as the proton cyclotron line, photo- ionization edge , 
and atomic tran sition lines at £ < 1 keV jHo et alj I2OO3I ; 
IPotekhin et ZI20041 However, hydrogen atmospheres cannot pro- 
duce strong spectral features at £ > 1 keV because (1) the 
binding energy of a hydrogen atom n ever exceeds ~ 1 keV at 
any B Q 1200 it ISanwal et al] l2002h . (2) QED effects signifi- 
cantly reduce line strengths at B > 7 x 10'^ G fao & Laiir2003l 
|2004| - Ivan Adelsberg & Laill200^ . and (3) the fraction of hydro- 
gen molecular ions (which may have transition lines above 1 keV; 
iTurbiner & Lopez Vi evra''2004^) is negligible (< 10"*) at B < 10'^ 
G and r > 10** K (Potekhin & Chabrier 2004). Therefore, spec- 
tral features at > 1 keV (such as the 1.4 keV absorption line of 
IE 1207) hint at a non-hydrogenic element atmosphere on the NS 
surfac e jSanwal et al1l2002l : iHailev & MorillioOa ; iMori & HailevI 
I2OO6I) . Also, recent theoretical studies suggest that a hydrogen layer 
on the surface may be rapidly depleted by diffuse nuclear burning 
or pul sar winds, thus exposing the heavier elemen ts that lie under- 
neath jChang et al.ll2004i ; IChang & Bildsteij26o4h . 

The few existing non-h ydrogenic atmosphere models l lMillerl 
I1992I : iRaiagopal et all [19971) are far from complete, due mainly 
to a lack of accurate atomic data for multi-electron ions in the 
high magnetic field regime and crude treatment of ionic motional 
and non-ideal effects. In this paper, we construct mid-Z element 
(carbon, oxygen and neon) atmosphere models that are compat- 
ible with the high quality X-ray data of Chandra and XMM- 
Newton. We only consider magnetic field strengths B = 10'"- 
10'^ G and efi"ective temperatures T^g = (1 - 5) x 10* K be- 
cause QED effec ts may become important at higher magnetic fields 
Jho & Lai|[2004l and references there in) and molecules m ay be- 
come abundant at lower temperatures l lMedin & Laill2006iJ) . QED 
effects on mid-Z elemen t atmosphere struc ture and spectra will be 
studied in future work l lHo&Morill2007h . We also defer inves- 
tigation of strongly-magnetized helium atmospheres due to non- 
trivial moti onal Stark effects an d possible contamination by helium 
molecules jMori & Hevill2007l) . 

We briefly describe here several major improvements of our 

mid-Z element atmosphere models over the previous atmosphere 
I — '. 1 (~~^ — I 

models by Miller (1992) (for helium and mid-Z elements) and 
[Raiagopal et al. ( 1997) (for iron). 

(1) The previous non-hydrogenic atmosphere models assume 
that bound electrons are in the ground Landau state (adiabatic ap- 
proximation). However, the adiabatic approximation is not valid for 
the low B and high atomic number Z regime that we consider in this 
paper. For instance, the ionization energy of the innermost electron 
of Z = 6-10 element ions at 5 = 10'^ G is underestimated by 
5-10% in the adiabatic approximation. 

(2) Numerous bound states and transitions are needed tocon- 
struct mid-Z element atmospheres. The model of iMilleJ i ll 9921) 
does not inclu de bound-bound transi tions, while the iron atmo- 
sphere model of lRaiagopal et al.l ( 1 19971) utilizes rather crude atomic 
line data (the energy values and oscillator strengths have as much 
as 10% and a factor of 2 uncertainties, respectively). We explicitly 
calculated line energies and oscillator strengths for numerous tran- 
sitions. For instance, nearly 600 transition lines are included in the 
oxygen atmosphere model at B = 10'- G. 

(3) In the case of hydrogen and helium, it is non-trivial to take 
into account the effects of a finite nuclear mass (Potekhin 1991; 
iBezchastnov et al .1 19981). However a simp le perturbative approach, 
introduced by IPavlov & MeszarosI jl993r) . is applicable to mid-Z 
element atmospheres due to the larger nuclear mass and binding 
energies. 

(4) All the previous magnetized atmosphere models, except 



the recent hydrogen atmosphere models of IPotekhin et al.l ( |2004| ). 
adopt polarization vectors assuming the plasma is fully-ionized. 
However, such atmosphere models may not produce correct spectra 
and absorption line profiles at low temperatures when the pla sma 
is partially-ionized jBulik & Pavlovll996l : |Potekhin et al.l2004l) . In 
our model, we explicitly compute polarization vectors by taking 
into account the effects of bound-bound and bound-free transitions. 

(5) Our radiative transfer code has been extensively tested for 
partially-i onized hydrogen atmospheres at vari ous magnetic field 
strengths i IHo et al.|[2003l ; IPotekhin et al]|2004l) . The code can be 
made applicable to mid-Z element atmospheres in which numerous 
line and edge features are present. We solve the coupled radiative 
transfer equations for the two photon polarizatio n normal mod es, 
unlike the polarization-averaged opacities used bv lMillej l fT993) . 

In this paper, we present each step in the construction of 
mid-Z element atmosphere models. We address important phys- 
ical effects in strongly-magnetized dense pl asmas, following the 

S iartial l y-ionized hydro gen calculations of Potekhin & Chabrieg 
20031) , IHo et alj l l2003h and IPotekhin et al] ( l2004h . Sections H 
and|4|discuss the atomic physics, equation of state, and opacities, 
respectiv ely, that we us e in th e m odelling; more details can be 
found in lMori& HailevI ( |2002|) and iMori & HailevI 1 I2OO6I) . Some 
aspects of the radiative transfer calculation are given in Section (H 
details of the radiative transfer scheme are described in IHo & Lail 
( I2OOII) and lHoetalJ l l2003h . Section|6|shows our atmosphere struc- 
ture results, and our model spectra are discussed in Section |7j we 
compare our mid-Z element atmosphere models with partially ion- 
ized hydrogen atmosphere models to illustrate the differences Q. A 
discussion on the determination of NS parameters is given in Sec- 
tionjSj We compare our model spectra to the observations of INSs 
in Section|9l SectionfTOlsummarizes our results. 



2 ATOMIC PHYSICS IN STRONG MAGNETIC FIELD 

In a strong magnetic field, the atomic structure is quite different 
from the non-magnetic case. A strong magnetic field deforms the 
atom into a cylindrical shape when magnetic field effects are larger 
than Coulomb field effects, i.e., at /Jz > 1 (Landau regime), where 
Pz = S/(4.7 X lO'Z^ G). A bound electron in the Landau regime 
is often denoted by two quantum numbers in and v [hereafter we 
use (mv) to denote a bound state], where in is the magnetic quan- 
tum number and v is the longitudinal quantum number along the 
field line, y = (tightly-bound) states have a larger binding en- 
ergy than y > (loosely-bound) states. When ySz 1> the ground 
state configuration is composed of bound electrons all in tightly- 
bound states, i.e. (mO) states. When the axisymmetry around the 
magnetic field is preserved, transitions with Am = and Ay = odd 
or Am = ±1 and Av = even are allowed by dipole selection rules. 
The former transitions occur in the longitudinal polarization mode 
(a = 0, parallel to 6), while the latter transitions occur in the cir- 
cular polarization modes (a = ±1, transverse to B). Several forbid- 
den transitions (e.g.. Am = 2 transition) will have non-zero oscilla- 
tor strengths when the moti onal Stark field breaks the a xisymmety 
around the magnetic field jPotekhin & Chabriej|2003h . However, 



' Most of the previous NS atmos phere models have been constructed by 
assuming a ful l y-ionized plasma jShibanov et al.l 1 19921 ; IZane et al] I2OO0I ; 
IHo & Laj|20"0ll ; IOzelll200l] ; lLlovj2003l) . However, the atomic fraction can 
be sizeable for hydrogen atmospheres at th e magnetic field strengths and 
tem peratures typicai of the observed INSs iHo et al J|2003| ; jPotekhin et alj 
I2QO4. 
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we neglect these weak transitions because the motional Stark ef- 
fects are significantly smaller in mid-Z element plasmas. 

We develop ed a fast and accurate atomic code suitable for the 
Landau regime jMori & Haiievll2002h . The code provides atomic 
data to better than 1% and 10% accuracy for energies and oscillator 
strengths, respectively. Most importantly, by applying perturbation 
methods to higher Landau levels, we have an extended range of 
applicable magnetic fields {fiz 0.3). The fast algorithm enables 
us to compute atomic data for numerous electron configurations 
within reasonable CPU time. We also evalu ate bound-free cross- 
sections following IPotekhin & Pavlov! ( 1 1997h . 



2.1 Line broadening 

Among various line broadening mechanisms in strongly- 
magnetized mid -Z element plasmas, motional Stark effects are the 
most important jMori & Hailevl200d) . As ions move randomly in a 
thermal plasma, coupling of the collective motion and the internal 
electronic structure induces a motional Stark field {Ems = 1^)5 
where ^ is the pseudomomentum and M is the mass of the 
ion. We evaluated the energy shift ca used by the motiona l Stark 
field using the perturb ation method jPavlov & Meszaroj 1 19931 ; 
ISezchastnov et all 19981) because the atomic mass is much heavier 
and the bin ding energies are significan tly larger than in the hydro- 
gen atom. IPavlov & Meszaro s[ ( Il993h introduced a perturbation 
method for hydrogen atom while the method of Bezchastnov et al] 
( 1 19981) is applicable to charged ions. In the perturbation method, 
the second-order energy shift (note that the first-order term 

vanishes) is given by 5^^' where A'j. is the pseudomomentum 
component perpendicular to the mag netic field and Mj^ ,. is the 
tran sverse mass (Mj^ ^. > M ) defined in lPavlov & MeszarosI i ll 9931) 
and lBezchastnov et alj ( Il998l) . 

For a given bound state k, one can evaluate Mx,i from atomic 
binding e nergies calculated under the infinite nuclear mass as- 
sumption jPavlov & Meszaroslll993l) . The increase in the mass due 
to motional Stark effects is roughly proportional to Z^^M"' for 
loosely-bound (v > 0) states. Therefore, the mass correction fac- 
tor of the H-like carbon ion is smaller than that of hydrogen by a 
factor of ~ 400. We computed the transverse mass for bound states 
with an electron in a v > state using the atomic data of H-like ions 
since bound electrons in a v > state are located at a large distance 
from the other bound electrons in tightly-bound states. This pro- 
cedure saves a significant amount of computation time, while the 
results are nearly identical to those of the more rigorous calcula- 
tion. Since the motional Stark field reduces th e binding energy, the 
line p rofile has a wing towards lower energies jPavlov & MeszarosI 
I1993I) . In our opacity calculation, we computed a line profile for 
mot ional Stark broadening fo llowing IPavlov & Meszaro3 jl993l) 
and lPavlov &Potekhinl ( ll995l) . 

Other less important sources of line broadening are thermal 
Doppler and pressure broadening. In a dense plasma, electron 
bound states are perturbed by electric fields from electrons (elec- 
tron collisional broad ening) and ions (quasi-static Stark broaden- 
ing) jSalzmannll 19981) . We evalu ated line widths by t he two types 
of pressure broadening following IPavlov & Potekhid l ll99^ . Pres- 
sure broadening becomes dominant only at high density (p > 
10^ g cm"^). 



2.2 Ionization balance 

In local thermodynamic equilibrium (LTE) atmospheres, the ion- 
ization fraction and level population are determined by the Saha- 
Boltzmann equilibrium. The validity of the LTE assumption will 
be discussed later (see Section |6^ . We iteratively solve the gen- 
eralized Saha-Boltzma nn equations in a magnetic field until con- 
vergence is achieve d jKhersonskiil 1 19871 : iRaiagopal et"aLl 1 19971 : 
iMori & Haiievll2006l) . We take into account pressure ionization by 
assigning occupation probab ilities for bound states using the elec- 
tric microfield distribution of lPotekhin et al.l 1 I2OO2I) . 

Figure [T] compares the eff^ective charge (Z^ff) of carbon, oxy- 
gen, and neon at 10'^ G and 10" G. The effective charge is defined 
as Zeff = where Z, and x, are the charge and the number 

fraction of ions of the ;-th type (2, x, = 1). We adopt the same tem- 
perature (r = 3 X 10* K) to illustrate the density dependence. In all 
cases, the degree of ionization decreases with density until pressure 
ionization becomes important at high density (e.g., p ~ 10^ g cm"^^ 
at fi = 10'^ G). Magnetized h ydrogen plasmas exhibit a similar be- 
havior jPotekhin etahll 19991) . though pressure ionization becomes 
important at lower densities in hydrogen atmospheres. 

The ion fractions have a strong dependence on magnetic field. 
The degree of ionization decreases at higher magnetic fields since 
the binding energy increases. Note that an unmagnetized carbon 
plasma is highly ion ized {Z^g ^ 5-6) even at low temperatures 
jPotekhin etai]|2005l) . At higher magnetic fields, pressure ioniza- 
tion takes place at higher density because the size of atoms and ions 
is smaller. 

At high magnetic fields a nd low tem peratures, molecules may 
be formed in the atmosphere ( lLaill200i]) . We find that neutral car- 
bon becomes abundant at T < 10' K and B = 10" G. Since we 
do not take into account any molecular transitions, our models be- 
come less accu rate in the regime whe re molecular contamination is 
non-neg ligible. iMedin & Lai ( l2006al) calculate binding energies of 
non-hydrogenic molecules using density-functional theory. Using 
their results, we solve the Saha equation fo r dissociation equilib - 
rium between the C atom and Ci molecule l lLai & SalDetei|[l997l) . 
Figure |2] shows the fractional ratio of the C atom and C2 molecule 
at 6 = 10" G. The solid curve is the density-temperature profile 
of a carbon atmosphere with B = 10" G and T^ff = 10* K (see 
Section |6j. Despite the fact that the internal degrees of freedom 
in molecules (such as ro-vibrational states) ca n enhance molecula r 
abundance by up to two orders of magnitude jMori & Hev'ill2007l) . 
we find that molecules can be safely ignored in our study. 



3 EQUATION OF STATE 

Given the results from Section |2j2| we compute the equation of 
state (EOS) of magnetized mid-Z element atmospheres with 



PI + K, + pz + PI + Pi- 



rn 



The first two terms {P'^ and P'.j) represent the ideal gas pressure 
of electrons and ions. Our calculation of P'^ takes into account 
the Landau quantization and electron degeneracy in a strong mag- 
netic field. In a dense plasma, non-ideal elfects are important, and 
they can induce negative excess pressure from electron-electron 
(P^j.), electron-ion (/"",.) and ion-ion coupling (P" , ). We compute 
the ex cess p ressure followingllchimaru et alJ ( ll987h , |Potekhin et aU 
( Il999l) . and lPotekhin & Chabrieil ( l2003h . 

Figure |3] shows the normalized pressure (P/NkT, where A' is 
the baryon number density) of a carbon plasma at different mag- 
netic field strengths and temperatures. We also plot the different 
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Figure 1. Effective charge Z^g as a function of density for carbon, oxygen, 
and neon plasmas at B = lO'^ G (top) and B = lO'^ G (bottom) and T = 
3 X 10* K. 

pressure components in Figure [S] The total pressure roughly fol- 
lows the ion fraction curves in Figure \T\ The total pressure de- 
creases with density until the plasma becomes more ionized by 
pressure ionization. At high densities (e.g., p > 10^ g cm"' for 
B = lO''^ G), electron degeneracy (dashed line) enhances the pres- 
sure sharply although it is partially compensated by the increasing 
excess pressure (dotted line). At higher magnetic fields, such tran- 
sition points shift to higher density because the binding energy is 
increased and the degree of electron degeneracy is decreased. 

For a given magnetic field strength and atmosphere compo- 
sition, we produced an EOS table in the temperature range 5.4 ^ 
log[o T ^ 7.8 with a step-size AlogjQ T = 0.3 and in the density 
range of -7.0 ^ logjoP ^ 5.0 with a step-size Alogmp = 1.0. The 
tables extend to low densities in order to cover strong absorption 
lines which form in the shallow layers of the atmosphere. 
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Figure 2. The fractional ratio of the C atom and C2 molecule at B = lO" G 
(dotted lines). The solid curve is the carbon atmosphere profile with B = 
10" G and T^fj = 10* K, as discussed in Section[6] 
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4 OPACITIES 

In this section, we briefly describe our opacity calculation. More 
details can be found in iMori & H ailev (2006). In a strongly- 
magnetized plasma, opacity becomes highly anisotropic, and ra- 
diation propagates in two normal modes: the ordinary mode (O- 
mode), which is mostly polarized parallel to the magnetic field, and 
the extraordinary mode (X-mod e), which i s mostly polarized per- 
pendicular to the magnetic field l iGinzburd[l97Q : iMeszaroSi,1992h . 



Figure 3. Normahzed pressure as a function of density for a carbon plasma 
at B = IOI2 G (top) and B = 10'' G (bottom) and T = 3 X 10* K. The solid, 
dashed and dotted are total, ideal (electron+ion) and excess pressure. Note 
that we plot the magnitude of the excess pressure since it is negative. 
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Figure 4. Cross-section (top), in units of Tliomson cross-section, and po- 
larizability tensor (bottom) as a function of energy in tlie a = mode for 
carbon at B = 10'^ G, T = 5 X lO' K andp = lO"^ g ^.^-i jjje dashed line 
shows the polarizability tensor for a fully-ionized plasma. The dominant 
ionization states are Li-like and Be-like carbon. 
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Figure 5. X-mode opacity of carbon at S = lO'^ G, T = 2 X 10* K, p = 
10"^ g cm"^, and 6 = 10°. The solid and dashed lines show the X-mode 
opacity for the complete and hybrid models, respectively. The fraction of 
bare and H-like carbon is 81 and 19%, respectively, while other ionization 
species have fractions less than 0.01%. The inset is a magnified view of the 
E = 0.4-0.7 keV region to illustrate the difference between the complete 
and hybiid models in the vicinity of the strongest bound-bound transition 
fine at ~ 0.6 keV. 



For a given density and temperature, the opacity (in cm- g ') in the 
normal mode / is given by 

K'{E,e) = Yj\<^i{E,efK„(E), (2) 

a 

where is the polarization eigenvector; |e/,| reflects the dielectric 
properties of the magnetized plasma and depends on the angle be- 
tween the magnetic field and the direction of photon propagation 6, 
such that = ei is the z-component (along B) of the mode eigen- 
vector and e{ = {e{ ± iei,)/ Vs are the circular components. 

Since mid-Z element atmospheres are usually not fully- 
ionized, we include the eff ects of bound species by using of 
the K ramers-Kronig relation JSulik & Pavlov|[l99d|Potekhin et al.l 
|2004) . Figure|4]shows the cross-section and polarizability tensor in 
the a = mode when carbon is partially-ionized (Zed = 2.9). The 
cross-section (ctq) is dominated by bound-bound and bound-free 
transitions in the X-ray band. The polarizability tensor (xo) devi- 
ates significantly from that in a fully-ionized plasma, particularly 
near the strong lines. 

When the degree of ionization increases, the difference be- 
tween using the polarizability tensors for partial ionization (com- 
plete model) and full ionization (hybrid model) becomes smaller. 
Figure |5] shows the X-mode opacity of carbon at B = 10'' G in 
the complete and hybrid models for a nearly fully-ionized case 
(Zeff = 5.8). Although the overall difference is tiny, there are no- 
ticeable deviations in the vicinity of the strong bound-bound tran- 
sition lines at E ~ 0.6 keV (see inset). Our results indicate that 
the strength and shape of strong absorption lines can be strongly 
affected even when the plasma is nearly fully-ionized. 

We produced opacity tables for the same temperature and den- 
sity grids as the EOS tables. For a given set of temperature and 
density, our opacity tables consist of cross-sections cr„(ii) and po- 
larizability tensors Xa(E) in the three cyclic polarization modes 
(a = 0, ±1). We adopted ~ 2500 equal logarithmically-spaced en- 
ergy grid points; this allows us to resolve narrow lines. 



5 RADIATIVE TRANSFER 

Thermal radiation from a NS is mediated by an atmosphere (with 
scale height ~ 1 cm) that covers the stellar surface. To determine 
the emission properties of a magnetized atmosphere, the radiative 
transfer equatio ns for the two coupled photon p olarization modes 
are solved (see IHo & Laill200ll ; IHo et al.ll2003l . for details on the 
construction of the atmosphere models). We note here a few partic- 
ulars of our models. 

The EOS tables for a given magnetic field are arranged by 
temperature and density. The pressure at Thomson depth t [= 
-K^' jpdz, where k^^ - {Z/A)cri/mp, the atomic mass number is 
A, and the Thomson cross-section is ctj] is calculated from hydro- 
static equilibrium. With a given temperature profile, we search the 
EOS tables for the nearest temperatures and gas pressures and do a 
weighted average to obtain the density profile. To include radiation 
pressure Pra d, we calculate P,-„ii from its contribution to hydrostatic 
equilibrium jMihalaslll978h 

dP,-ad o-sBr'ff 
(It c k^^ ' 

where 

/ = i^X ^'^^ 

Fj is the specific flux for mode j, and F is total flux; this contri- 
bution is usually small (see Section IgTI . The radiation pressure is 
then subtracted from the total pressure to obtain the gas pressure. 

The self-consistency of the atmosphere model is deter- 
mined by requiring that the fractional temperature corrections 
AT(t)/T(t) <g: 1% at each Thomson depth, deviations from ra- 
diative equilibrium « 1%, and deviations from constant total flux 
< 1%. Note that the atmosphere models formally have a depen- 
dence, through hydrostatic balance, on the surface gravity g [= 
(1 + Zg)GM/R^] and thus the NS mass M and radius R; however, 
the emergent spectra do not vary much using differ ent values of 
g around 2 x lO'"* cm s"^ (see Section 1731 see also IPavlov et al.l 
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Figure 6. Temperature profiles of carbon (solid lines) and liydrogen atmo- 
spheres (dashed lines) with B = 10'^ G and Teff = 3 X 10* K. The dotted 
lines are grey profiles, r is the Thomson depth and tr is the Rosseland 
depth. At large r, temperature profiles merge with grey profiles. 

1 19951) . Most of our models are constructed using a surface grav- 
ity g = 2.4 X lO"* cm s"- with M = lAM^, R = W km, and 
Zg = 0.3 (see Section l731 for other cases). Also, though our atmo- 
sphere models can have a magnetic field at an arbitrary angle ©g 
relative to the surface normal, most of the models considered in 
Section|7]have the magnetic field aligned perpendicular to the stel- 
lar surface (0g = 0°; see Sections [7.3| and [F!2l for other cases). The 
spectra represent emission from a local patch of the NS surface. 



6 ATMOSPHERE PROFILE 

Figures [6] and |7] show the temperature profiles of carbon atmo- 
spheres in comparison with those of partially ionized hydrogen at- 
mospheres at 6 = 10'^ G and B = 10'^ G. There are two distinct 
diflrerences between hydrogen and mid-Z element atmosphere pro- 
files; these difi^erences are due to the fact that bound-bound and 
bound-free opacities are more important than free-free opacities. 
Temperature profiles in the mid-Z element atmospheres are closer 
to grey than in hydrogen atmospheres since the opacities are less 
energy-dependent, while the free-free opacity (which is dominant 
in hydrogen atmospheres) has a steep energy dependence. Photo- 
sphere densities are also significantly smaller since opacities are 
higher than in hydrogen atmospheres, and photosphere densities 
increase with magnetic field since opacities are reduced. 

Figure [8] shows the effective charge of carbon atmospheres at 
different optical depths. At B = 10'^ G and T^f{>3x 10« K, the H- 
like carbon ion is largely populated. At higher B and lower T^if, var- 
ious ionization states exist in the atmosphere, and they are likely to 
produce numerous spectral features. Figure [9] illustrates where X- 
mode and O-mode photons decouple from the matter as a function 
of energy. It is clear that X-mode photons come from deeper, hotter 
layers in the atmosphere than O-mode photons; therefore spectra 
are strongly polarized and dominated by the X-mode photons. 
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Figure 7. Temperature profiles of carbon and hydrogen with B = lO'^ G 
and Teff = 3 X lO'' K. Notation is the same as in Figure|6] 
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Figure 8. Elfective charge of carbon as a function of Thomson depth r at 
B = 10'2 G and B = lO'^ G and T^g = 10^ and 3 X 10* K. Carbon is nearly 
fully-ionized at Teff = 5 X 10* K. 

6.1 Radiative levitation 

At high temperature, mid-Z element ions may be levitated by ra- 
diation pressure through strong absorption lines. Our scheme de- 
scribed in Section [5] does not perform radiative transfer of indi- 
vidual lines since we compute the gradient of radiation pressure 
using opacities that are integrated over frequency. Although de- 
tailed line-by-line radiative transfer modeling is beyond the scope 
of this paper, we estimate here the effects of radiative levitation in 
mid-Z element atmospheres. Gravitational force on a carbon ion 
is 4.6 X 10"' dyne when g = 2.4 x lO'"* cm s"^, while radiation 
force through an absorption line is 1.6xlO""(T/(rjy//10' K)** dyne, 
where ct/ is the line cross-section in units of Thomson cross-section 
(cTr). In the carbon atmosphere at fi = 10'^ G and Jeff = 5 x 10* K, 
several strong lines have cross-sections as large as ~ 10"* 0-7-. There- 
fore, radiative levitation can be important through strong absorption 
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Figure 9. Thomson depth (where photons decouple from matter) as a func- 
tion of energy for a carbon atmosphere with B = 10'- G and T^ff = 
3 X 10^ K. The bottom two panels plot the local temperature and density 
at the coiTesponding Thomson depth. The solid and dotted lines are for 
the X-mode and O-mode, respectively. For the two polarization modes, the 
photon-decoupling layer is defined such that the effec tive (angle - averag ed) 
optical depth of each mode is unity. See Section 4.2 in lHo & Laa )200lh for 
more details on how we calculate effective optical depths. 
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lines at T^g ~ 5 x 10 K, while it will be negligible at lower tem- 
peratures. 

6.2 Validity of LTE assumption 

The LTE assumption may become invalid at low density and high 
temperature. To our knowledge, the validity of LTE assumption 
has not been investigated in previous atmosphere models for cool- 
ing NSs. Three conditions need to be fulfilled for a LTE plasma 
zmann|[l998l) : (1) Maxwell-Boltzmann distribution of the elec- 
tron and ion velocities, (2) Saha equilibrium for the charge state 
distribution, and (3) Boltzmann distribution of the excited states. 
The first condition is well-satisfied in our case since the electron 
self-collision time is very short. However, the latter two conditions 
are usu ally more restrict ive. We examine the second condition fol- 
lowing ISalzmanij l fl993) . For the third condition, we estimate ra- 
diative and coUisional transition rates between the g round and ex- 
cited states and apply the criteria from lGrienil ( Il964l) . 

In Figure[TO] we plot the (p, T) regime where the LTE assump- 
tion is valid, along with carbon atmosphere profiles at 6 = 10'^ G 
and two effective temperatures. The condition for a Boltzmann dis- 
tribution (hatched area) is more restrictive than the one for Saha 
equilibrium (grey area). We find that both conditions break down in 
the upper atmosphere layers with p < 10"* g cm"' and 10"' g cm"' 
at T^g = 10* K and 5 x 10* K, respectively. Therefore, only a 
few of the strongest lines will be affected by non-LTE effects. At 
higher densities, Saha equilibrium is recovered. The condition for 
the Boltzmann distribution (hatched area) becomes irrelevant in 
the (p, T) regime where a fraction of excited states is small (solid 
lines). Therefore, the LTE assumption is completely valid above 



Figure 10. LTE conditions for carbon atmospheres at B = lO'^ G and 
Teff = 10* K (top) and T^g = 5x 10* K (bottom). The Saha equation is valid 
in the grey area, while the Boltzmann distribution for the excited states is 
valid in the hatched area. The solid + dashed curves show the temperature 
profile of the atmospheres. The fraction of excited states exceeds 10% on 
the dashed lines, while the occupation probability of the ground state is 
larger than 90% on the solid lines. 



p ~ 10"* g cm"' at Teff = 10* K. On the other hand, the occupation 
probabilities of excited states computed by using the Boltzmann 
distribution may not be correct in the region between p ~ 4 x 10"' 
and 3 X 10"^ g cm"' at Teff = 5 x 10* K. We find that such a region 
does not exist at lower effective temperatures. 

We thus expect that only a few strong absorption lines will be 
altered by non-LTE effects, while the bulk of the mid-Z element at- 
mosphere structure and spectra are unaffected in our models. To 
address the non-LTE effects, one must find collisional-radiative 
steady states by solving rate equations with b oth radiative and col- 
lisional processes included l lSalzmanrill 19981) . Such an analysis is 
beyond the scope of this paper. 



7 MODEL ATMOSPHERE SPECTRA 

We construct spectra of mid-Z element atmospheres with B = 
10'" G and 10" G and three different effective temperatures (Teff = 
10*, 3 X 10*, and 5 x 10* K). The results are shown for carbon (Fig- 
ures [TT] and [T2]l, oxygen (Figures [T3l and [T4t and neon (Figure [TS] 
andll6t. For comparison, we plot partially ionized hydrogen atmo- 
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Figure 11. Spectra of carbon atmospheres with B = 10 G and three ef- 
fective temperatures (T^ff = 1,3 and 5 X 10^ K). The dashed lines show 
partially ionized hydrogen atmosphere spectra with the same B and T^g, 
and the dotted lines show blackbody spectra with T = T^g. 
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Figure 12. Spectra of carbon atmospheres with B = 10 G. Notation is the 
same as in Fig. II II 

sphere spectra (dashed lines) and blackbody spectra (dotted lines) 
with the same magnetic field strengths and temperatures. 

7.1 Spectral energy distribution 

Overall, mid-Z element atmosphere spectra are significantly softer 
than hydrogen atmosphere spectra (dashed lines) and close to 
blackbody spectra (dotted lines) because temperature profiles are 
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Figure 13. Spectra of oxygen atmospheres with B = 10 G. Notation is 
the same as in Fig. II II 
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Figure 14. Spectra of oxygen atmospheres with B = 10 G. Notation is 
the same as in Fig. II II 

closer to grey. The mid-Z element atmosphere spectra are even 
softer than blackbody in the Wien tail due to the photo-absorption 
edges (at ~ 1.5-2 keV for B = 10'^ G models). The absorption 
edges are due to bound-free transitions from the innermost elec- 
tron [(mv) = (00) state] in different ionization species. At higher 
energy, the spectra become harder than blackbody and hydrogen 
atmospheres in order to maintain radiative equilibrium. A simi- 
lar feature is seen in the magnetized iron atmosphere models of 
iRaiagopal et al.l l ll997l) . This appears to be a common signature of 
heavy element atmospheres over a large range of B and T^g. 
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Figure 15. Spectra of neon atmospheres with B = lO'" G. Notation is the 
same as in Fig. 1111 
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Figure 16. Spectra of neon atmospheres with B = W'^ G. Notation is the 
same as in Fig. 1111 



7.2 Absorption features 

Numerous absorption lines are present, especially in low tempera- 
ture models, and as expected, heavier element atmospheres show 
more absorption features. Magnetic field dependence is also re- 
markable. There are more discrete absorption lines at B = 10'^ G 
than at S = 10'^ G. Most of the absorption lines at fi = 10'^ G 
are longitudinal transitions (Am = 0, Av > in the a = mode). 
On the other hand, the londitudinal transitions become significantly 
weaker at higher B because the a = component of the polarization 



vector components in the X-mode becomes smaller with increas- 
ing B. This effect is most noticeable in the low temperature models 
(e.g., see the neon models at T^ff = 10* K). The suppression of 
longitud inal transitions is also seen in magnetized hydrogen atmo- 
spheres bio et alj|2003h . Note that longitudinal transitions appear 
strong in the mid-Z element atmospheres at fi = 10'^ G since the 
magnetic field is not as dominant over the nuclear Coulomb field at 
j8z> 1. 

In the B = 10'^ G model spectra, it is mainly the transverse 
transition lines that are visible. At S = 10''' G, Am = -l-l,Av = 
transitions from bound electrons in outer shells are strong (e.g., 
lines at ~ 0.26, 0.39, 0.69 and 1.8 keV in oxygen atmosphere mod- 
els at i? = 10'^^ G). The lines appear at lower energies when the 
temperature is lower and oxygen is less ionized (compare oxygen 
spectra at T^g = K and 3 x 10* K for B = 10'' G in Figure[T4Tl. 
Other weaker and broader transverse transition lines {Am = ±1, 
Ay > 0) are also present in the B = 10" G model spectra. 

Most absorption lines have asymmetric profiles with low en- 
ergy tails due to motional Stark effects. Some of the strong lines 
show more complex profiles due to the enhanced polarizability near 
the lines. In general, line widths become larger at higher magnetic 
fields since motional Stark effects are more important. Only photo- 
absorption edges from the innermost electrons are clearly visible at 
high energy where there are no bound-bound transition lines (e.g., 
£ > 1.7 keV for carbon atmosphere spectra at 6 = 10'^ G), al- 
though bound-free transitions are important throughout the X-ray 
band. 

7.3 Dependence on magnetic field geometry and surface 
gravity 

Here we illustrate the dependence of the atmosphere models on 
magnetic field geometry (&b) and surface gravity (g). Figure [TT] 
shows temperature profiles and carbon atmosphere spectra at two 
different magnetic field angles (Ob) with respect to the NS sur- 
face normal. There is no remarkable difference between the two 
magnetic field orientations, except at several strong lines and the 
photo-ionization edges at high energy (E > 1.7 keV for carbon at 
B = 10'- G). The angular dependence shown here reflects the (rota- 
tional) phase variation of atmosphere spectra, as the observer sees 
regions of the NS surface with different magnetic field orientations 
(see Section ing . 

Atmosphere emissi on can also be highly beamed in the pres - 
ence of a magnetic field jShibanov et alj|l992l : |Pavlov et al.lll994l) . 
Figure [TH] shows the energy-dependent beam pattern for the at- 
mosphere models with T^ff = 3 x 10* K and B = 10'^ and 
10" G. Note that the emission is more beamed at higher magnetic 
fields. There are some noticable differences between hydrogen and 
mid-Z element atmospheres. Our results show that energy-resolved 
lightcurve analysis can potentially be used as a diagnostic tool to 
determine the surface composition. 

Figure [19] shows carbon spectra at two rather extreme values 
of g (= 5 X 10 '^ and 5 x lO'"' cm s~^) an d hence NS mass and ra- 
dius (see, e.g.. lLattimer & Prakashll200lh . Our results show mid-Z 
ele ment atmosphere sp ectra, like magnetic hydrogen spectra (see, 
e.g.. |Pavlov et alj|l995l) . have a weak dependence on gravitational 
acceleration. This is contrary to the results of unmagnetized NS at- 
mosphere spectra in which several a bsorption features v ary in their 
strengths and shapes strongly with g jZavlin et al.ll996[) . The weak 
dependence on g is partially due to the fact that spectral features in 
mid-Z element atmospheres are primarily broadened by motional 
Stark effects (which are independent of plasma density) and pres- 
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Figure 17. Temperature profile (top) and carbon atmosphere spectra (bot- 
tom) with &B = 0° (solid lines; magnetic field perpendicular to the at- 
mosphere plane) and &b = 90° (dashed lines; magnetic field parallel to 
the atmosphere plane). The other parameters are set to B = lO'^ G and 
Teff = 3 X 10« K. 



sure effects are small (since absorption lines are formed at low den- 
sity). 



8 X-RAY SPECTROSCOPY AND DETERMINATION OF 
NEUTRON STAR PARAMETERS 

Given the results in the previous sections, it is useful to consider 
what one can learn from X-ray spectroscopy of INSs, especially 
those which show absorption features in their spectra. Four NS 
parameters are important in determining thermal spectra; surface 
composition (Z), magnetic field (6), effective temperature (Tcb), 
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Figure 18. Polar diagrams of the specific intensity (scaled to the intensity 
in the normal direction) from the carbon (solid lines) and hydrogen (dashed 
lines) atmosphere models with T^g = 3 X 10* K and B = 10'^ G (top) and 
B = 10'3 G (bottom). 



and gravitational redshift (Zg). It is desirable to detect multiple spec- 
tral features to break degeneracies between the parameters. We dis- 
cuss here the feasibility of measuring these NS parameters with 
X-ray spectroscopy. 



8.1 Surface composition 

Given the substantial differences in spectral energy distribution 
(SED) and absorption features, it should be relatively easy to distin- 
guish between mid-Z element and hydrogen atmospheres. In partic- 
ular, a possible signature for the presence of a mid-Z element atmo- 
sphere is the detection of multiple strong absorption lines at high 
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Figure 19. Temperature profiles (top) and cai'bon atmosphere spectra (bot- 



tom) with 



5 x 10" (dotted fines), 2.4 X 10" (solid fines), and 



5 X lO''' cm s ^ (dashed lines), corresponding to {M,R) = (0.5Mo,12 km), 
(1.4Mo,10 km), and (2.2Mo,10 km), respectively. The other parameters are 
set toB= 10'2 G, Teff = 3 X 10'^ K, and = 0°. 



energy (E > I keV). Multi- wavelength SED analysis over optical 
and X-ray bands will be useful since mid-Z element atmospheres 
are si gnificantly softer than hydrogen atmospheres jPavlov et al.l 
Il996h . Since the number of absorption features and their energies 
are strongly dependent on the surface composition, it is possible 
to distinguish between carbon and oxygen atmospheres if multiple 
spectral features are detected. 
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Figure 20. Carbon atmosphere spectra with B = 2 X lO'- G and &b = 0° 
(solid line) and B = lO'- G and ©g = 90° (dashed line). In both cases, 
Teff = 3 X 10'' K. The former and latter cases correspond to atmosphere 
spectra from the pole and equator in a dipole magnetic field configuration, 
respectively. 



8.2 Magnetic field 

The absorption features shown in Figures [TTt|16| are d eeper and nar- 
rower than the observed spectral features from INSs ("Sanw al et al.l 
12002; Haberl et al. 2003, 2004b; van Kerkwijk et al. 2004). This is 
partly due to the fact that our atmosphere models are constructed 
with a single magnetic field strength B and orientation 0g, thus re- 
flecting emission from a local patch on the NS surface. The absorp- 
tion features in our model spectra will appear broadened if there 
are magnetic field variations over the surface. We expect to see 
dramatic changes in the spectra even with a small magnetic field 
difference, since the ionization fraction and line energies depend 
strongly on B. Figure |20] shows carbon atmosphere spectra at two 
magnetic field strengths and angles, corresponding to the pole and 
equator for a magnetic dipole configuration. Indeed, atmosphere 
spectra can look quite dilferent merely due to a factor of two differ- 
ence in the magnetic field strength. Thus thermal spectra of INSs 
can show significant rotational phase-variation due to a magnetic 
field distribution, even when the surface temperature is uniform. 

In order to construct synthetic atmosphere spectra, one must 
integrate emission from patches of the NS surface with differ- 
ing B, @B, and r^ff. Such synthetic spectra are necess arily model- 
dependent (see, e.g.,[Zavlin et al. 1995; Z ane et al.l200 1; Ho & La^ 
l2004IZane & Turollalbooel ; IHo et al.ll2007l) . as the magnetic field 
and temperature distributions over the surface are unknown. Al- 
though magnetic field broadening is different for each absorption 
line, atomic line energies vary approximately by AE/E ^ 10% 
and 30% for a factor of 1.5 and 2 variation in magnetic field 
strength. Electron or ion cyclotron lines may appear in the X-ray 
band when B is lower or higher than the values considered in 
this paper. Some of the obs erved absorption fe atures have been 
attributed to cyclotron lines jHaberl et al.ll2003l) . Since cyclotron 
line energy varies linearly with B, we expect cyclotron lines to be 
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broader and have larger phase variation than atomic lines. There- 
fore, line width analysis can be a useful diagnostic tool to probe 
the magnetic field distribution over the NS surface, given that the 
intrinsic widt hs of absorption lines are narrow in mid-Z element 
atmospheres l lMori & Haiievll2006h . Magnetic field variation can 
smear out narrow absorption lines and reduce the accuracy of de- 
termining other NS parameters such as the gravitational redshift. 
Phase-resolved spectroscopy is essential to disentangle magnetic 
field effects from the other parameters and to probe the magnetic 
field geometry. 



8.3 Effective temperature 

The effective temper ature must be measured to constrain the cool - 
ing history of NSs jPage et al.ll2004l : lYakovlev & Pethickll2004t) . 
Thermal spectra of most INSs are fit with either blackbody or fully- 
ionized hydrogen atmosphere models, each of which exhibits the 
softest and hardest spectra, respectively. Therefore, errors in the 
measured effective temperatures are primarily due to uncertainty in 
the surface composition. 

Although mid-Z element atmospheres can have overall SEDs 
close to a blackbody, there are several significant differences from 
blackbody spectra. In the Wien tail, mid-Z element atmosphere 
spectra are softer than a blackbody due to photo-absorption lines 
and edges and become harder than a blackbody at higher energy. 
Therefore the effective temperature can be higher than the black- 
body temperature if high energy tails are not fully covered by X- 
ray spectroscopy. Alternatively, the high energy part of X-ray spec- 
tra may be erroneously fit by an additional blackbody component 
(see Section |9T| for the case of 1E1207). It is important to cover 
a wide energy band and fit the data with model atmosphere spec- 
tra, rather than fitting with one or two blackbody components and 
simple Gaussian or Lorentzian absorption lines. 



8.4 Neutron star mass and radius 

One of the ultimate goals in NS research is to uniquely determine 
the mass and radius of NSs and constrain the nuclear equation of 
state of the NS core. In X-ray spectroscopy, gravitational redshift 
Zg is one of the NS parameters which can be measured once X-ray 
data are fit with model atmosphere spectra. The accuracy of a grav- 
itational redshift measurement is limited by magnetic field surface 
variations since line energies change with both B and Zg- To in- 
crease accuracy, phase-resolved spectroscopy is essential because B 
varies with rotation phase but Zg does not. Grating spectroscopy will 
measure gravitational redshift with high accuracy if magnetic field 
variation is not so severe and photon statistics is good. In the fu- 
ture, high-resolution spectroscopy missions, such as Constellation- 
X, will be able to measure the magnetic field and gravitational red- 
shift with unprecedented accuracy. 

It is possible to measure the surface gravity (g) from line 
widths if the lines are primarily broadened by pressure effects 
jPaerelj Il997h . However, motional Stark effects produce line 
broadening far larger than pressure effects; this is partly because 
lines are formed in low density layers. Also, our study indicates that 
it is difficult to measure surface gravity from overall spectral fitting 
since mid-Z element spectra have a very weak dependence on g 
(see Figure [T9t. Therefore, the methods for constraining g by using 
absorption features is applicable only to weakly-magnetized NSs 
(B < 10' G) since their spectra are nearly fre e from the magneti c 
field effects discussed in this paper (see also IZavlin et al. 



On the other hand, it may be possible to constrain g in heavier ele- 
ment atmospheres (e.g., Fe) since motional Stark effects are further 
reduced while ion-coupling effects are enhanced. We plan to inves- 
tigate heavy element atmospheres composed of silicon and iron in 
future work. 

Radius measurements are feasible by combining X-ray spec- 
troscopy (which measures flux) and distance measurements. Dis- 
tances to several nearby INSs have been m easured with rela- 
tively high accuracy by parallax measurements jWalter & Lattimed 
l200ilKaDlan et al.l2002l : lvan Kerkwiik & Kaplanl20"o^ . However, 
it should be kept in mind that only a fraction of the surface (such as 
hot spots) may be visibl e due to temperature anisotropy caused by 
strong mag netic fields jGreenstein & Hartkel 1 1 9831 : iGeppert et al.l 
I2004ll200a) . Independent lightcurve analysis will be helpful to de- 
termine the size of the emission region. 



9 COMPARISON WITH X-RAY DATA OF ISOLATED 
NEUTRON STARS 

In this section, we discuss INS thermal spectra in the context of 
mid-Z element atmospheres. We consider the INS 1E1207.4-5209 
and the class of radio-quiet INSs; these all show spectral features in 
their X-ray spectra, except RX J1856. 5-3754. Table[T]summarizes 
their observational properties based on recent Chandra and XMM- 
Newton observations. 



9.1 1E1207.4-5209 

1E1207 is a hot isolated with an age ~ 7x lO"* y. Remarkably, it 
shows two broad absorption features at ~ 0.7 and ~ 1.4 keV and an 
excess above blackbody emission at > 2 keV. A second blackbody 
component with a small emitting area (i? < 1 km) and high temper- 
ature (r > 3 x 10^ K) has been used to fit the high energy excess 
jMereghetti et alj|2002l ; ISignami et alll2003l ; lOe Luca et al]|2004 
iMori et al.l2005l) . However, it is puzzling that the pulsed fraction is 
not particularly large at high energies (De Luca et al. 2004), while 
small hot spots can strongly modulate X-ray lightcurves. Also, the 
temperature can be quite different depending on the number of con- 
tinuum components used to fit the data, e.g., kTBB=260 eV when 
using a single blackbody and kTBB=l50 eV when using a double 
blackbody (the s econd blackbody in the latter case has ~300 eV; 
iMori et al.ll2005l) . 

Proposed models involving cyclotron li nes or atomic lines 
from a light element (H, He) atmosphere jSanwal et alj 1 20021 ; 
' r^ ' 



iBignami et al.l2003l ; lTurbiner & Lopez Vievral2004l ) seem unlikely 
since line str engths are significantly weakened by vacuu m reso- 
nance effects jHo & Laill2004l ; lvan Adelsberg & Lail2006l. and ref- 
erences therein) or the ioniz ation stat es responsible for the ob- 
served lines are not ab und ant I Mori & H ailev 200^. Alternatively, 
lHailev&Moril ( l2002l) and lMori & HailevI ( i2006i) proposed a mid-Z 
element atmosphere for IE 1207. 

Figure [2T| shows the XMM-NewtonfEFlC-'P'N spectrum of 
1E1207 and the oxygen atmosphere spectrum with B = 10'- G, 
Teff = 4 X 10^ K and Zg = 0.4; the latter is convolved with the 
XMM-Newton EPIC energy resolution. Although many of the nar- 
row absorption lines are blended or smeared out, some substructure 



^ Recent timing a nalysis suggests that IE 1207 is in a binaiy system 
iWoods et al]|2006h . However, it is unlikely that 1E1207 is an accreting 
NS. 
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Object 


kTBB [eV] 


Bd [1013 


£•,,„, [keV] 


(T [eV] 


EW [eV] 


1E1207.4-5209*''' 


260 


7 


0.7, 1.4 


130, 76 


100, 60 


RBS 1774 


102 


< 24 


0.75 


27 


27 


RX J1308.6+2127 


103 


3.4 


0.23, 0.46 


150, 260 


200, 180 


RXJ1605.3+3249 


96 


7 


0.40, 0.59, 0.78 


87'" 


96, 76, 67 


RXJ0806.4-4123<''> 


96 


< 14 


0.43 


70(f) 


33 


RXJ0720.4-3125'"' 


85-95 


2.4 


0.28 


90 


0-70 


RX J1856.5-3754 


62 


^ iW) 








RX J0420.0-5022<=' 


45 


< 18 


0.33 


70(f) 


45 



Notes: 

We list blackbody temperature (Tbb), dipole magnetic field strength (Bj), line energy, line width, and equivalent width of the observed absorption features. All 
the temperatures were measured by assuming a sing le blackbody comp onent. The observed absorption lines were fit with a Gaussian line p rofile and cr refers 
to the Gaussian line width de fined in equati on (1 ) oflMori et al.1 j2005h . The line parameters are subject to the choice of continuum models iMori et ani2005h 
and the number of fitted lines iHaberl200d) . See lHabe'm200^ for a complete reference list of the observations. 

(a) The source shows an excess above a blackbody at £ > 2 keV; an additional continuum component has been used to fit this excess iMori et alj2005h . 

(b) An additional absorption feature may be present at higher energy, and spectral fitting with two lines yields £"/,;,(. = 0. 3 1 and 0.6 1 keV iHaberj|2006l). 

(c) Over the last six years. X-ray spectra of RXJ0720 have become harder, and the line strength has increased ide Vries et al J l2004t I VinkTrakl |2004| : 
iHaberl et alj2006l) . .— ,, ,— , 

(d) Dipole magnetic field strength was estimated from t he Ho' nebula in the v icinity of RX Jl 856.5-3754 Ivan Kerkwiik & Kulkaml200il : lKaDlan et alj2002h . 

(e) Based on the 20 ksec XMM-Newton EPIC-PN data, lHaberl eTall j2004al) indicated an absorption feature at £ = 0.33 keV. However, its significance is not 
definitive yet (Saberi 2006), and our independent analysis of the XMM-Newton data has not confirmed a significant spectral feature. 

(f) The line width was fixed while fitting the spectra. 



is still seen. Remarkably, the oxygen atmosphere spectrum in Fig- 
ure [21] exhibits the three spectral features seen in the 1E1207 spec- 
trum (i.e., two broad absorption features at ~ 0.7 a nd 1 .4 keV and 
spectr al hardening above ~ 2 keV). See Table 1 in lMori & HailevI 
( l2006l) for identification of the bound-bound transitions responsible 
for the 0.7 and 1.4 keV features and their unredshifted line ener- 
gies. In our convolved model spectrum, the 0.7 keV feature shows 
a more complex structure than the 1 .4 keV feature because there are 
more blended absorption lines at ~ 0.7 keV. Indeed, some analyses 
of the XMM-NewtonlEPlC data of IE1207 indicate the presence 
of sub struture at ~ 0.7 keV jMereghetti et alJl2002l : lMori & HailevI 
l2006h . The edge-like feature in the model spectrum at ~ 2 keV 
(similar to the excess above a blackbody in the IE 1207 spectrum) 
is due to bound-free transitions from the innermost electron in the 
(mv) = (00) state of various oxygen ions. In the 1E1207 spec- 
tra, there is an indication of an additional absorption feature at 
E < 0.3 keV, where the oxygen atmosphere spectrum also shows 
an absorption line. Given the less reliable instrumental calibration 
and stronger effects of interstellar absorption at low energies, we 
plan to more carefully investigate the significance of this potential 
absorption feature in future work. 



10-" r 



Note that the model atmosphere spectrum shown here is not 
completely realistic because we have not taken into account the 
distribution of magnetic field and temperature over the surface (see 
Sections [7 . 3 1 and [F!2l ) . More detailed analysis based on spectral fit- 
ting by varying NS atmosphere parameters and viewing geometries 
is in progress. Give n the good photon statistics with the 2 60 ksec 
XMM-Newton data l lBignami et alJl2003l : lM"ori et alj|2005l) and the 
three (or possibly four) spectral features to fit, we are optimistic 
about pinning down the surface composition, magnetic field, tem- 
perature, and gravitational redshift of 1E1207 with the existing 
CCD spectral data. 
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Figure 21. Unfolded XMM-Newton EPIC-PN spectrum of 1E1207.4-5209 
(top) and spectra of oxygen atmosphere with B = 10'^ G, Teff = 4 X 10* K 
and Zg = 0.4 (bottom). In the top panel, we overlaid a blackbody spectrum 
to illustrate the absorption features; we also show the ratio of the data and 
blackbody. There is potentially another absorption feature at £ < 0.3 keV 
(indicated by an arrow in the ratio plot). In the bottom panel, the solid line 
is the model atmosphere spectrum convolved with the XMM-Newton EPIC- 
PN energy resolution, while the dashed line is the raw spectrum. 
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9.2 Radio-quiet neutron stars 

Seven nearby INS are characterized as radio-quiet NSs (RQNSs) 
since tliey exhibit similar spectral and timing properties 
( Ivan Kerkwiik & KaplanI l2006h . At least five of them show a 
single or multiple absorption features at £ 0.2 - 0.7 keV 



jHaberl et alj|2003l; Ivan Kerkwiik et al 


l2004lHaberl et alj|2004hl: 


IZane et alj 20051: ISchwooe et alj 2005 


: Haberll 2006h. while the 



brightest INS RX J1856.5-3754 shows a r iearly p erfect black- 
body spectrum without any spectral features jBurwitz et al. 20o3)- 
Some of the observed absorption features have been interpreted 
as proton cy clotron lines or hydrogen/heliu m atomic lines at 6 = 
lO'^-IO''* G ( I van Kerkwiik & KaDlanl2006|). Timing analyses yield 



similar dipole magnetic field strengths i Kaplan & van KerkwiikI 
l2005al lbl). 

Figure [22] shows the XMM-NewtonfEFlC-'P'N spectrum of 
RX JI605. 3+3249 and the oxygen atmosphere spectrum with B = 
10" G, Teff = 1.5 X 10^ K and Zg = 0.4; the latter is convolved 
with XMM-Newton EPIC energy resolution. We use the spectrum 
of RX JI605. 3+3249 as an illustration of the X-ray spectra of sev- 
eral other RQNSs with absorption features at £ ~ 0.3-0.4 keV (see 
Table [TJ. In addition to the feature in the model spectrum at (red- 
shifted) 0.4 keV, there are weaker absorption features at higher en- 
ergies; these may be pr esent in the sp ectrum of RX JI605. 3+3249 
or RX JI308.6+2127 l lHaberll 120061) . If all RQNSs have similar 
magnetic field strengths (e.g., B =^ few x lO'^ G, as suggested by 
timing analyses), we may expect RQNSs with higher temperatures 
to have more distinct absorption features at high energy (see, e.g., 
Figure[T6]l. This appears to agree with the observed properties listed 
in Table [T] Further X-ray observations are needed to confirm the 
significance of multiple absorption lines. 

On the other hand, it is puzzling as to why the spectrum of 
RX J1856.5-3754 can be fit by a featureless blackbody spectrum. 
Since RX JI856. 5-3754 has a lower temperature than most of the 
other RQN Ss, it is possible that molecular chains are formed on 
the surface jMedin & Laill2006alB), or RX J 1856.5-3754 has the 



condensed surface i Turolla et al.fl2004l ; Ivan Adelsberg et al]|2005l ; 



IPerez-Azorm et al.ll2005l) . Recent theoretical studies indicate that 
the condensed surfaces of magnetized NSs exhibit blackbody-like 
spect r a in the X-ray band jTurolla et alJI 20041 ; Ivan Adelsberg et al] 
I2OO5I ; IPerez-Azorm et al.ll2005l) . At the observed temperature of 
RX J1856.5-3754 {Tbb ~ 7 x 10^^ K), it is estimated that a car- 
bon adno^spher£jmay_becoj^ when 6 > 3 X 10'^ G 
( lLaill200ll : [Medin & Lail2006bl ). In this context, we expect that the 
other RQNS (RX J0420.0-5022) with a low temperature {Tbb ~ 
5 X 10"* K) has a condensed surface; RX J0420-5022 should thus 
show a featureless blackbody-like spectrum U unless it has a sur- 
face magnetic field lower than 10'^ G. Indeed, the sharp tempera- 
ture dependence of the molecular abundance was shown in the case 
of a helium atmosphere jMori & Hevl2007l) . 

Optical flux excess from Rayleigh-Jeans tails of X-ray 
blackbody spectra is another common signature of RQNSs 
teurwitz et alj [2OO3I) . Fully-ionized hydrogen atmosphere mod- 
els overpredict the optical fluxes by 1-2 orders of magnitude 
( |Pavlovetal.lll996l ; | Pons et al] 1 20021) . Unmagnetized heavy el- 
ement atmosphe re models can account for the optical excess 
jPons et alj[2002h . but they also exhibit numerous absorption fea- 
tures, which significantly deviate from the featureless X-ray spec- 
tra of RX JI856.5-3754 jZavlin et al.lll996l ; lRaiagoDal & Rom.mil 



^ See caption (e) of table 1 for the current status of the putative line feature 
in the RX J0420.0-5022 spectrum. 
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Figure 22. Unfolded XMM-Newton EPIC-PN spectmrn of 
RXJ1605. 3+3249 (top) and spectra of oxygen atmosphere with B = iO'^ G, 
Tejf = 1.5 X 10* K and z.g = 0.4 (bottom). In the top panel, we overlaid a 
blackbody spectrum to illustrate the absorption features; we also show the 
ratio of the data and blackbody. In the bottom panel, the solid line is the 
model atmosphere spectrum convolved with the XMM-Newton EPIC-PN 
energy resolution, while the dashed line is the raw spectrum. 



ll996l ; lGansicke et al.ll2002h. Theref ore, it requires either a thin hy- 
drogen atmosphere 1 Ho et al.ll2007h or an anisotropic temperature 
distribution over the surface l|Ponsetal.i r2002: Braie & RomanJ 
l2002l ; |pirez-Azorfn et al.ll2006l) to produce both the featureless X- 
ray spectrum and the optical excess. Since magnetized mid-Z ele- 
ment atmospheres have SEDs close to blackbody, they may account 
for the optical excess. Unlike unmagnetized atmospheres, magnetic 
field variations can smear out absorption features in magnetized 
mid-Z atmosphere spectra and produce featureless X-ray spectra. 



10 SUMMARY 

We have constructed equation of state and opacity tables for mid-Z 
elements using the most up-to-date methods available and improv- 
ing upon earlier calculations by including more physical effects. 
Using these tables, we built magnetic neutron star atmosphere mod- 
els for carbon, oxygen, and neon. Mid-Z element atmosphere spec- 
tra show some dependence on magnetic field orientation (relative to 
the surface) and a weak dependence on surface gravity. There are 
many characteristics that distinguish these spectra from hydrogen 
atmosphere spectra, e.g., a multitude of lines, different line energy 
and strength dependences on the magnetic field, and possible hard- 
ening at high energies. 

There are still several issues that need to be addressed in order 
to improve our models. These include accounting for molecules. 
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the quantum electrodynamical effect of vacuum resonance, and 
non-LTE effects. Molecules and vacuum resonance may be impor- 
tant for atmosphere models with magnetic fields higher than those 
considered in this paper, i.e., B » 10'^ G JHo & Morill2007h . We 
also plan to study atmosphere models with different elements, such 
as helium, silicon and iron, in order to cover a larger NS atmo- 
sphere parameter space. We will construct synthetic atmosphere 
model spectra by taking into account magnetic field and temper- 
ature distributions on the surface. 

We presented qualitative comparisons between our model 
spectra and the observations of several isolated neutron stars. The 
results are promising; the model atmosphere spectra closely resem- 
ble the X-ray spectra of these objects. The similarities strongly 
motivate more detailed studies, including phase-resolved analysis, 
since line identification and a better temperature estimation can 
provide invaluable information on neutron stars. 
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